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Abstract 



I We show that with a few modifications the Adomian's method for 

I solving second order differential equations can be used to obtain the 

' known results of the special functions of mathematical physics. The 

, modifications are necessary in order to take corectly into account the 

behavior of the solutions in the neighborhood of the singular points. 



1 Introduction 



Recently a great deal of interest has been focused on the application of Ado- 
mian's decomposition method to solve a wide variety of stochastic and de- 
^ ' terministic problems [^.Although the Adomian's goal is to find a method 

■ to unify linear and nonlinear, ordinary or partial differential equations for 

solving initial and boundary value problems, we shall deal in the following 
only with linear second order differential equations. 

Our aim here is to compare the decomposition method with the classical 
methods for solving differential equations in order to obtain a better under- 
standing of it. Because the special functions are extremely useful tools for 
obtaining closed form as well as series solutions to a variety of problems aris- 
ing in science and engineering we tryed to reobtain the known results by the 
new method. We found that with a slight modification the method works 
also in this case. 

Usually one starts with an equation Ff{x) = g{x), where F represents 
a general nonlinear ordinary differential operator. F is decomposed in two 
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parts L + R, where L is the hnear part of F . The core of the method 
consists in finding an inverse of L . In many cases this is not possible, in 
others the Green function is quite comphcate and consequently there are 
difficulties of integration. The solution proposed by Adomian is to take 
L as the heighest-ordered derivative of the linear part. For example, for 
initial value problems Adomian defines for L = (P /dx^ as the two-fold 
integration operator from to x . We have L~^Lf{x) = f{x) — /(O) — x/'(0) 
and therefore 

fix) = /(O) + xfiO) - L-'R + L-'g 1.1 

Such stated the method is quite old. For linear second order differential 
equations it has been proposed by Cochran, who transformed the equation 

\p{x)f'{x)]' + q{x)f{x) = g{x) 

into an equivalent Volterra integral equation . Cochran has been interested 
in the uniqueness of solutions for two-point boundary-value problems and in 
the existence of eigenfunctions of the related homogeneous equation. 

Adomian's point of view was to take L as the simplest easily invertible 
operator pushing all the other terms in the reminder R or in the non- homo- 
geneous part. With this choice the method has given many beautiful results. 

Our first comment here is that in order to find the classical results for the 
special functions this approach is not sufficient. The second one is that we 
give an example, the Bessel polynomials equation, where L has to be chosen 
as a first order differential operator for solving the equation, although the 
equation satisfied by them is of second order. The last comment points out 
that the approach 1.1 supposes that f{x) and its derivative are finite at the 
point X = , thus one cannot obtain solutions in the neighborhood of a 
singular point. 

As a first conclusion this means that the decomposition method is a useful 
tool for solving concrete problems, but clearly its form cannot be decided 
independent of the problem. 

We shall consider the linear second order equation written in the more 
general form 

L{x, D)f{x) - R{x, /, Df, D'f) = 1.2 

where D = d/dx and R is the reminder which includes all the other terms 
not contained in the principal part L{x, D) whose form is 

L{x,D) = h{x)Dp{x)D 1.3 

h{x) and p{x) are smooth functions whose properties are given below. The 
minus sign in 1.2 is taken for convenience. 
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At this stage the clue of the problem is to make a clever choice of L and 
it! such that the resulting pseudo Volterra integral equation should be easily 
solvable. 

In many cases bringing 1.2 to the canonical form complicates the matter. 
Wc propose to find an inverse for the general form 1.3 and use it to find 
many of the known formulas for the solutions of classical special function 
equations. 

Our second approach consists in writting 1.2 as a formal non-homogeneous 
equation and by using a variant of the method of variation of parameters we 
transform it into a pseudo Volterra integral equation. We consider this form 
as the most general form of the decomposition method for the second order 
differential equations. We apply it to solve the Bessel functions equation. 



2 Special Functions Equations 

Most of the equations for the classical special functions are of the form 1.2, 
with the principal part L{x, D) of the form 1.3. Many of them have h{x) = 1. 

The functions h{x) andp(a;) are supposed to be smooth and l/p{x) locally 
integrable around some point, which we can take to be x = without loss of 
generality. The decomposition method consists in finding an inverse operator 
and with its help we are able to construct in many cases the exact solution 
of the problem. 

A formal inverse of 1.3 can be easily found. We choose it as 

One sees that 

{LL-')f{x) = I-f{x) 

but 

iL-'L)fix) ^ I ■ fix) 

where / is the identity operator. The second relation tell us that L^^ is not a 
true inverse, it becomes so when we take into account the initial conditions, 
which is equivalent with prescribing the values of /(O) and f{0) . 
Indeed 

dt 



fix)-m-pio)f'{o) [ 

Jo 



Pit) 
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Thus we obtain from 1.2 

m = m - p(o)/'(o) I'^^i'^j; ^fiv, m 2.2 

which is a Vohcrra integral equation.. Its most general form for second order 
differential equations will be given in Section 3. The solution of the last 
equation is sought as a series 

00 

k=0 

where 

fo^f{o)-p{o)m 



Pit) 

and fk{x) is obtained by the Picard method of succesive approximation 

fk^L-\x,D)fk-i 

In the following we shall apply the above technique to a few classical 
equations. 

2.1 Gauss Hypergeometric Function 

Because the solutions to many mathematical physics problems are expressed 
by the hypergeometric functions we consider first this equation whose stan- 
dard form is 

x{l - x)f"{x) + [7 - (a + /3 + l)x]f{x) - a(3f{x) = 
and rewrite it as 

Lf{x) = xf"{x) + jf'ix) = x'-^^{x^^) 

= x'^f'ix) + {a + (3 + l)xf(x) + a(3f(x) 
In this case h{x) — x^~'^ and p{x) — x'^. We suppose that 7 > and then 

{L-'f){x)=f{x)-f{0) 

We take as usual in the decomposition method 

f{x)=foix)+Mx)+f2{x) + ... 

and choose the initial condition as 

/o = /(O) = 1 
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Using the form 2.2 we get the recurrence relation 

x~''dx / x'^~^[x'^ f^_i{x) + {a + P + l)xf^_i{x) + a(3fn-i{x)]dx 

JO 

We find that 

aP X 

Jl = TT 

7 1! 

and by induction 

. _ (a)n(/5)nX" 
In 



(7)n 

where {a)n = «(« + 1) . . . (a + n — 1) 

In this way we get the well known solution 

fix) = ,F,ia, /3; 7; X) = 1 + £ i^^^ 

n=l W/n ^• 

2.2 Confluent Hyper geometric Functions 

As concerns the confiuent hypergeometric equation 

xf'ix) + (7 - x)f'{x) - af{x) = 

its solution can be obtained in the same manner as that of the Gauss hyper- 
geometric function, the functions h{x) and p{x) being the same. 
We shall consider now the equation 

f"{x) + ax^-^f'ix) + bx^-^fix) = 2.3 

where a, b and q are complex numbers, whose solution seems to be unknown. 
For g = it reduces to Euler equation. For others values of q it is equivalent 
with a degenerate hypergeometric function. More about this equation see P, 
eq. 2.60]. 

We construct two independent solutions in the neighborhood of the point 
X = 0. We take L = D"^ and then 

L-V(x) = /(x)-/(0)-x/'(0) 

The first solution is defined by 

/(O) = 1, /'(O) = 2.4 

We choose 

fo = /(O) = 1 
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and find 



We make the ansatz 

fk = Ckx'"^ 

and from tlie relation f^^i = L~^fk we get 



JO Jo 

_ aqk + b x^''+^^i 

~ ~'"'{k + l)q-lq{k+l) 
In this way we find the recurrence relation 

aqk + b 

'^^^ = ~[q{k+l)-l]q{k+lf' 
Taking into account 2.5 we find 

^ ^ (g-l)(2g-l)...(ng-l) 



^ ^ (l-l/g)„U/ 
The solution of the Eq. 2.3 with the initial conditions 2.4 is 

f(x) = i + f(-ir (-Y— 

i.e. f{x) is the confluent hypergeometric function 

/ x = 2.6 

ag g g 

The second solution which satisfies 

/(O) = 0, /'(O) = 1 
is found in the same way and is given by 

J., . ^ ,a + 6 g + 1 ax^, 

f{x) = X iFi , ^ ; 2.7 

aq q q 

The solutions 2.6 and 2.7 give the general solution of Eq.2.3 and as func- 
tion of q have poles at g = ±l/n + 1, respectively, which accumulate at g = 0. 
When h/aq = —p, and a + b/aq = —n, respectively the solutions reduce to 
polynomials in x*. 
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2.3 Bessel Polynomials 

Until now in all the applications of the decomposition method the operator 
L contained a piece involving the highest-ordered derivative. We want to 
present an example where such an approach does not work. If we form L 
from lower-ordered derivatives we can solve the equation. 
We consider the equation 

x^f\x) + {2ax + P)f\x) - 7(7 + = 2.8 

where a, 13 and 7 are complex numbers. If a = 1, /? = 2 and 7 = n the 
above equation is the defining equation for the Bessel polynomials 0. 
We take L as the first order operator 

dx 

The above equation is rewritten as 

Lf{x) = fix) = i[7(7 + I) fix) - 2xaf'ix) - x'f'ix)] 

In this way we find the Volterra equation 

fix) = /(O) + ^ r[7(7 + I) fix) - 2xaf'ix) - x'fix)] dx 
p Jo 

We look for a regular solution of Eq. 2.8 at the origin and take the initial 
condition 

/o = /(O) = 1 



We get 



We make the ansatz 



^ 7(7 + 1) X 
p 1! 



^ Cfc 



From the integral equation we find 

^ Ck[e + (2a - l)k - 7(7 + 1)] x^+^ 
ik + l)\ 

Now we define 

(2z/+ 1)2 = (2a -1)2 + 47(7 + 1) 
and from the above relations one finds the recurrence relation 



Cfc+i = -(/c + a + z/)(/c + a - V - l)ck 
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whose solution is 

Cfc = i-l)^ia + z/)fc(a -iy-l)k 
and the solution of Eq. 2.7 is 

X 

f{x) = 2Fo{a -u -l,a + u; --) 
For a — V — 1 = —n we find the generalized Bessel polynomials [4]. 

3 Volterra Integral Equation 

Let be a general second-order differential equation 

f{x) + a{x)f'{x) + b{x)f{x) = h{x) 3.1 

and let he ip{x) ^ a. solution of the corresponding homogeneous equation. 

By applying a variant of the method of variation of parameters we obtain 
the general solution of Eq. 3.1 as 

where E{x) = exp J a{x)dx and Ci , C2 are constants 0. 

In the above form is defined as an indefinite integral; for every prob- 
lem we have to transform it into a definite integral according with the solution 
we are looking for. 

We consider 3.2 as the starting point in the formulation of the most 
general form of the decomposition method. We use Eq. 3.1 to transform a 
linear homogeneous equation into a pseudo Volterra integral equation. The 
crucial point is to separate a convenient part in the left hand side of 3.1 such 
that we can easily find a solution of this part; this separation has to take into 
account the behavior of the solution in the neighborhood of the point where 
we are looking for the solution, to improve our chance to obtain the full series 
in explicit form. The simplest separation is by moving the term b{x)f{x) on 
the right-hand side, then a solution of the left side alone is /(x) = 1 . We 
shall give a few examples rather than treat the general case. 

Let us come back to the hypergeometric equation 2.1 and write it as 

fix) + ^f'ix) = xf'ix) + {a + (3 + l)nx) + ^/(x) 

X X 

We treat the right-hand side of this equation as a formal non- homogeneous 
term and look for a solution of the homogeneous part, i.e. of the left part. 
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The solutions are 93+ (x) = 1 and 93- (x) ~ x^~~^ We consider first </?+ and 
take Ci = 1 and C2 = . Since we look for the regular solution at x = we 
take as a 2-fold integration from to x and find from 3.2 the Volterra 
integral equation 

fix) = 1 + / ^ / x-^ixfix) + {a + /3+ l)f'{x) + ^f{x)]dx 3.3 
Jo x^ Jo X 

We solve Eq. 3.3 by iteration taking the zero-order approximation the non- 
homogeneous term /(O) = /o = 1 and find 

fX fjq- fX (yfi 

/n+1 = / -7 / ^"[^/:(^) + {a + f3 + + ^Ux)]dx 

J X J X 

which coincides with the result found in the preceding section. 

If we start with the second solution (/?_ and take as above Ci — 1 and 
C2 = we find the integral equation 



fX fX (yR 

f(x) = x'~^ + x'"^ / dx x''-^ / [xf(x) + (a + p+ l)f(x) + —f(x)]dx 
Jo Jo X 

We solve it by iteration, but now the zero-order approximation is /o = x^~'^, 
and find 

f _ i-, (a-7 + l)(/3-7 + l) X 
^'"^ 2-7 1! 

By induction we get the series for the second solution 

/(x) =x^-^2i^i(«-7 + l,/3-7 + l;2-7;^) 
We consider now the equation of the Bessel functions and write it as 

fix) + -fix) + (1 - -^)fix) = 

X X'^ 

Its solution cannot be obtained by the decomposition method as originally 
stated, nor in the more general form of it given by us in Sec. 2. The reason 
is the singular behavior of the Bessel functions in the neighborhood of the 
point X = 0; however its solution can be easily obtained using the form 3.2 . 
Indeed, we rewrite the Bessel equation in the form 

fix) + -fix) - -Jix) = -fix) 

X X 

and observe that the "homogeneous" part has the solutions (p±ix) — x^''' . 
Since the behaviour of the Bessel functions at x = is that given by (fi± we 
choose as a 2-fold integration from to x . 
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Let us try the first solution; we choose Ci = 2 ^ /Viv + 1) and 6*2 = 
and from 3.2 we get the integral equation 

\.{V + I) JO JO 

The initial approximation of the solution is the non-homogeneous term which 
is no more constant 

{x/2r 

r(i/+i) 

and from the above equation we get 

(.x/2)-+^ 

■'^ r{u + 2) 

and by induction we find the known series. 

fix) = Mx) = ^ ^ ^ 



+k+i) 

If we choose the initial conditions as 

2-u 



Ci = 0, C2 = 



r(i - u) 



proceeding as above we find the second independent solution, J_^(x) . 

The equation 3.2 represents the starting point for obtaining the most 
general form of the decomposition method. Indeed it can be recovered, in 
the form given by Adamian, by the choice a{x) = b{x) = in Eq. 3.1. As a 
consequence ^p{x) = 1 and E{x) = 1 and if we look for regular solutions in 
the neighborhood of a; = 3.2 has the form 

f[x) = Ci + C2X + dx h(x)dx 
Jo Jo 

and by identification 

Ci = f{0) C2 = f{0) 



4 Conclusion 

We have shown that with the necessary modificatons the Adomian's method 
can be used to obtain the classical results on the special functions. Rather 
than prescribe a unique form we show that the concrete problem decides 
how to formulate the method. The clue of this one consists in transforming 
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a linear differential equation into a pseudo Volterra integral equation whose 
solution is obtained by the Picard process of succesive approximation. 

The method can be obtained by a light modification of the known form of 
the solution for a linear second order differential equation, the sole difference 
being the extension of what is ussualy called the non-homogeneous term in 
the sense that it can include important pieces from the homogeneous part. 

Our results can also be seen as a good illustration for the effectivness of 
the Picard method of succesive approximation. 
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